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We present a theoretical study of an escape rate for switching from the superconducting state 
to a resistive one in series arrays of strongly interacting Josephson junctions. At low temperatures 
such a switching is determined by macroscopic quantum tunneling (MQT) of a spatially extended 
Josephson phase. An increase of the crossover temperature from the thermal to quantum regimes, 
and a giant enhancement of MQT escape rate is obtained. Such an effect is explained by excitation 
of collective charge oscillations forming a charge instanton. Using a model of screened Coulomb 
interaction we found that our analysis is in a good accord with recently published experimental 
results on an enhancement of the MQT in single crystals of high-T c superconductors. 
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Great attention has been devoted to an experimen- 
tal and theoretical study of diverse macroscopic quan- 
tum phenomena, e.g. macroscopic quantum tunneling 
(MQT) and energy level quantization just to name a few, 
in Josephson coupled systems 0, EI1I1I1I10I10. 
A lumped Josephson junction is characterized by a sin- 
gle Josephson phase, and the MQT of a Josephson phase 
has been found long time ago in Nb Josephson junction 
0. At low temperatures the MQT determines the 
escape rate of the switching from the superconducting 
state to a resistive one. The MQT is considered as a 
basic effect in the modern field of quantum information 
processing 01 

Next step in a study of macroscopic quantum phenom- 
ena is to obtain the MQT in spatially extended Joseph- 
son systems, i.e. dc SQUIDs |3j, quasi-one-dimensional 
long Josephson junctions 0, UTS 01' Josephson paral- 
lel and series arrays 0, . These systems that can be 
characterized by coordinate and time dependent Joseph- 
son phase ip(x,t), present a particular case of interacting 
many particle systems. The presence of inductive inter- 
action in Josephson parallel arrays and long Josephson 
junctions allows one to form diverse macroscopic objects 
(solitons), e.g. Josephson vortices (magnetic fluxons) and 
vortex-antivortex pairs. Thus, the MQT in spatially ex- 
tended Josephson systems can be often reduced to the 
quantum fluctuation induced escape of solitons from a 
pinning potential 0111101. Indeed, both the MQT of 
a Josephson vortex ||j and the quantum dissociation of 
a vortex-antivortex pair have been observed. 

Other example of spatially extended Josephson sys- 
tems is a dc biased series array of Josephson junctions. 
This case presents a special interest because artificial se- 
ries arrays of Al / AI2Q3 / Al junctions have been prepared 
and layered high-T c superconductors can be mod- 
elled as a stack of intrinsic Josephson junctions [l^ . 
Moreover, modern fabrication technique allows to pre- 
pare single crystals of layered high-T c superconductors 
with an extremely homogeneous distribution of critical 



currents of intrinsic Josephson junctions, and a low level 
of dissipation 

00. 

Such systems are suitable for an ex- 
perimental study of macroscopic quantum phenomena, 
and indeed, the MQT of a Josephson phase has been ob- 
served in layered BiSrCaCuO superconductors 00. 

A rather unexpected result found in the Refs. 0, 
is that the crossover temperature T* from the thermal 
fluctuation regime to the MQT regime is much larger in 
respect to a single Josephson junction having the same 
parameters. Moreover, in Ref. it has been observed 
that the MQT escape rate Tmqt for a stack of intrin- 
sic Josephson junctions is jour order of magnitude larger 
than r mqt for a single Josephson junction. In these ex- 
periments it was also found that the escape rate Tt in 
the thermal fluctuation regime did not differ from the 
escape rate of a single Josephson junction. 

Since in the model of independent Josephson junctions 
the crossover temperature does not depend on number 
A^ of junctions, and the escape rate T is just propor- 
tional to N, an enhancement of the MQT observed in 
layered high-T c superconductors stems from an interac- 
tion between intrinsic Josephson junctions. In order to 
explain such a giant increase of the MQT escape rate the 
idea of a long-range coupling between intrinsic Joseph- 
son junctions has been proposed in Ref. 0. Although 
this mechanism can definitely take place, we notice that 
a natural source of interaction in series junction arrays 
is the screened Coulomb charge interaction [l3|. This 
interaction is especially strong in layered high-T c super- 
conductors, where the Debye screening length is of the 
order of superconducting layer thickness. 

In this Letter we will show that the excitation of collec- 
tive charge oscillations in a stack of Josephson junctions 
leads to an increase of the crossover temperature T* and 
a giant enhancement of the MQT escape rate. In or- 
der to quantitatively analyze the escape rate in layered 
superconductors a series array of interacting Josephson 
junctions biased by dc current / (see schematic in Fig. 1) 
is considered. The system is characterized by the set of 
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FIG. 1: Schematic of a dc biased layered high-T — c supercon- 
ductor and a series array of Josephson junctions. Collective 
charge oscillations (wave lines), a strongly localized instan- 
ton (dashed line) and a charge instanton with long tails (solid 
lines) are shown. 



Josephson phases Lp n (r), where number n changes from 
to N. In order to obtain the escape rate T we use "in- 
stanton technique" [lil ITol Il6| , and therefore, r is the 
imaginary time that varies from to fijk^T (T is the 
temperature). In the imaginary time representation the 
Lagrangian of a series array with a charge interaction is 
written as 

L = i[v3„(r)C„ m ( / 9 m (r)] +^f7„ , 



U n {ip) = cos<p„(r) + itp n {r) , 



(1) 



where i = I / I c is the normalized external dc current, and 
I c is the nominal value of the critical current of a single 
junction. The Lagrangian is expressed in units of Ej, 
where Ej is the Josephson energy of a single junction. 
The capacitance matrix C nm = C|„_ m | is determined by 
a charge interaction between junctions, and C nm becomes 
the S nm - function for the specific case of decoupled junc- 
tions, i.e. as the charge interaction between junctions is 
small C << Co- The Co is the capacitance of a single 
Josephson junction. In the experiments on the MQT in 
Josephson coupled systems Ej >> Ec = e 2 /C*o, and 
the MQT occurs as the dc current I is close to I c , and 
therefore (i-l)< 1. In this case the potential U n ((p) is 
written as 



U n (ip) = (1-£)<^„(t)- 
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(2) 



The escape rate is determined by the particular solution 
fn (j) providing the extremum of action 

, r h/k B T 

S = ~hJ L{T)dT ■ (3) 

At high temperatures such a solution is determined by 
extremum points of the potential U n , and it is written as 

<p n = 2^2(1 - i)S nl - y/2(l-i) . 



Here, I is a junction number where the fluctuation occurs. 
Since this solution does not depend on the time r, we 
can immediately conclude that the excitation of charge 
oscillations does not change a slope of the dependence of 
the escape rate ln[IY(l — i)] on the bias current (1 — i) in 
the thermal fluctuation regime. However, the crossover 
temperature from the thermal fluctuation regime to the 
MQT can be strongly enhanced by excitation of charge 
oscillations in a stack of interacting Josephson junctions. 
Indeed, using the method elaborated in 0, 0, 0] we 
obtain that at high temperatures the optimal fluctuation 
(y9„(r) around an extremum point has a form: 



(4) 



where the eigenfunctions <p n are the solution of the non- 
local and inhomogeneous equation: 



{\-Lol)^ n . (5) 



Here, A are the eigenvalues of the Eq. 0, ojq = u p [2(l — 
i)] 1 / 4 is the dc bias dependent frequency of oscillations 
on the bottom of potential well, U n (ip), and lu p is the 
plasma frequency of a single Josephson junction. The 
crossover temperature T* is determined by the condition 
that there is the eigenvalue A = 0, [3]. Using the 
Fourier transform of Eq. (pj the crossover temperature 
is obtained as a solution of the following equation: 



1 = 
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2n Wq 



£ 2 [l + 2C(g)] 



2nk B T* 



(6) 



Here, C(q) is the Fourier cos-transform of the function 
C nm , i.e. C(q) = -^Y, n =iCnCOs(qn). We can ana- 
lyze Eq. © in different limits. In the case of decou- 
pled junctions, i.e as C(q) << 1 for all q, we obtain 
that T* nd — hu>o/ (2irkB), and the crossover temperature 
of a series array T* nd coincides with the one for a sin- 
gle Josephson junction. However, for strongly coupled 
Josephson junctions C(q) becomes negative for q ~ ir, 
and can be of the order of one. In this case, the crossover 
temperature has to increase. 

This effect can be easily understood if we notice that 
the spectrum of charge oscillations in a stack of Joseph- 
son junctions forms a band 



u 2 (q) = 



1 + 2C(q) 



(7) 



The crossover temperature is determined by the maxi- 
mum value of the oj{q) . Therefore, the excitation of high- 
frequency charge oscillations leads to an increase of the 
crossover temperature. On other hand a well-known ef- 
fect of an enhancement of the MQT escape rate under 
microwave radiation is determined in Josephson junc- 
tion series arrays with a strong charge interaction by the 



3 



minimum value of the uj q . Therefore, the value of the 
plasma frequency found from this effect is not directly 
related to the crossover temperature. 

Now we turn to the MQT regime, where the extremum 
point of the action S is the " tau-dependent" instanton 
(bounce) solution. In the absence of a charge interaction 
the instanton solution is strongly localized on a particular 
junction (see schematic in Fig. 1, dashed line), i.e 



Vnir) = /o(r) = S n l 



3^/2(1 -i) 
cosh 2 (loqt/2) 



(8) 



In a generic case of interacting junctions , an instanton 
dressed by collective charge oscillations has small tau- 
dependent spatial tails (see Fig. 1, solid line). Such a 
spatial-temporal instanton solution that can be called 
a charge instanton, is found by perturbation analysis 
(C(q) < 1) as 



f n {t) 



2k 



dq 
2^ 



dujdri 2C{q)t 



j(t-ti) 



2tt LU 2 +L0 2 [l + 2C(q)] 



/o(ri) 



(9) 

Substituting J^J in the expression for action S we 
obtain the escape rate Tmqt (in physical units) as 



Tmqt = r exp 



where 



X = 60-7T 



271 



dq 
2^ 



T2Ej_ 



dx- 



2 V4(l_ i) 6/4 (l _ x) 



C 2 (q)x 6 



, (io) 



{.T 2 [l + 2C(g)] + l}sinh^(7rx) 



(11) 

The parameter \ having a positive value, characterizes an 
enhancement of MQT due to the presence of the charge 
interaction between Josephson junctions, as C{q) is dif- 
ferent from zero. The integral over x can be calculated 
by taking into account that the typical values of x are 
numerically small, i.e x < (1/tt). Therefore, with a good 
accuracy we obtain 



X 
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y 



2tt 



C 2 (q) 



dq 
2^ 



(12) 



Thus, remarkably the crossover from the thermal fluc- 
tuation regime to the MQT regime in Josephson series 
arrays with a charge interaction is realized as a transi- 
tion from 0- to 2-dimensional instanton solution. 

Presented above a generic analysis can be used for a 
particular case of the MQT in a single crystal of high-T c 
superconductors contained intrinsic Josephson junctions. 
A natural source of a charge interaction in layered super- 
conductors is the Coulomb screened charge interaction, 
and we use a simple model introduced by M. Tachiki et 
al. in Ref. In this model the function C(q) has a 

form: 



C(q) 



a 



1 — cos(g) 



2 1 + a(l - cos(g)) ' 



(13) 




FIG. 2: The dependencies of the crossover temperature T*(a) 
and the parameter x( a ) on the strength of Coulomb interac- 
tion a. 



where the parameter a characterizes a strength of the 
charge interaction. This parameter is determined by the 
ratio between the Debye screening length and a supercon- 
ducting layer thickness, and a can be larger than one in 
layered high-T c superconductors [13J. Substituting 
in 10 and (|12l) . and calculating the integrals over q, we 
obtain the dependence of the crossover temperature T* 
and the parameter \ on the strength of Coulomb inter- 
action a: 



T*(a) = T t 



and 



X 



ind 




l + 3a 
(1 + 2a) 3 / 2 



(14) 



(15) 



These dependencies are shown in Fig. 2. 

Moreover, in Fig. 3 we present the dependence of the 
MQT escape rate on the external dc bias for various val- 
ues of a. One can see a giant enhancement of the MQT 
escape rate as the parameter a becomes larger than one. 
This enhancement results from a decrease of the slope 
of the bias current dependence escape rate ln[r^(l — i)] 
on (1 — i). Comparing our theoretical predictions with 
the experimental curves published in Ref. (see Fig. 5 
in Ref. 0) we find a good agreement as a ~ 4. By 
making use of the same value of a we obtain that the 
crossover temperature increases two times in respect to 
a single Josephson junction with the same parameters. 

Notice here that for a finite system as a number of 
intrinsic Josephson junctions N is not very large, there 
is a discrete number of collective oscillation modes that 
can be excited in a system. In this case the integral 
over q in Eq. (|12|l has to be changed to the sum over 
discrete modes, i. e. J dq/(2n) — > (l/iV)^ n . There- 
fore, the slope of the bias current dependence escape rate 
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FIG. 3: The dependence of the MQT escape rate Tmqt on the 
dc bias current 7 for various values of a = 0, 1, 3, 4. Legend 
corresponds to the curves arranged from right to left. In order 
to fit our theoretical predictions to the experimental results 
of Ref. the value of Ej/(Tiujp) ~ 264 and preexponent in 
the expression 11101 To = 1 GHz were used. 



ln[lY(l — i)] increases for a series array with less number 
of Josephson junctions. It gives a natural explanation of 
a moderate enhancement of the MQT escape rate as a 
number of intrinsic Josephson junction TV was decreased 
□ 

In conclusion we have shown that the MQT escape rate 
can be drastically increased in strongly coupled Joseph- 
son junctions series arrays. The MQT in such arrays is 
described through the quantum fluctuation induced exci- 
tation of a charge instanton characterized by long spatial 
tails. It is in marked contrast to the well-known cases of 
a single Josephson junction or weakly coupled Josephson 
junctions series arrays where a strongly localized instan- 
ton solution forms. Appearance of the charge instan- 
ton is explained by excitation of collective charge oscilla- 
tions, and this mechanism has to be especially effective 
in layered high-T c superconductors where the screening 
of Coulomb interaction is rather weak. The excitation 
of high-frequency charge oscillations also leads to a sig- 
nificant increase of the crossover temperature from the 
thermal fluctuation to the quantum regime. Such an en- 
hancement of the MQT in Josephson series arrays with 
a strong charge interaction can be very promising for a 



field of quantum information processing . 

I would like to thank Yu. Koval, P.Miiller, A. Volkov 
and A. V. Ustinov for useful discussions and attraction 
of my attention to the MQT problem in high-T c super- 
conductors. I acknowledge the financial support by SFB 
491. 



[1 
[2 

[3. 

[4 
[5 

[6] 

[7] 

[8 

[9 
[10 

[11 

[12 

[13 

[14 

[15 
[16 



M. Tinkham, Introduction to Superconductivity, 2nd edi- 
tion, McGraw-Hill, New York (1996) 
R. F. Voss and R. A. Webb, Phys. Rev. Lett. 47, 265 
(981); J. M. Martinis and J. Clarke, Phys. Rev. Lett. 55, 
1908 (1985). 

F. Balestro, J. Claudon, J. P. Pekola, and O. Buisson, 
Phys. Rev. Lett. 91, 158301 (2003); S.-X. Li, Y. Yu, Y. 
Zhang, W. Qiu, S. Han, and Z. Wang, Phys. Rev. Lett. 
89, 098301 (2002) 

T. Kato and M. Imada, J. Phys. Soc. Jpn. 65, 2963 
(1996); T. Kato, J. Phys. Soc. Jpn., 69 2735 (2000). 
A. Schnirman, E. Ben- Jacob, and B. A. Malomed, Phys. 
Rev. B. 56, 14677 (1997). 

A. Wallraff, A. Lukashenko, J. Lisenfeld, A. Kemp, Yu. 
Koval, M. V. Fistul, and A. V. Ustinov, Nature 425, 155 
(2003). 

M. V. Fistul, A. Wallraff, Yu. Koval, A. Lukashenko, B. 
A. Malomed, and A. V. Ustinov, Phys. Rev. Lett. 91, 
257004 (2003). 

K. Inomata, S. Sato, Koji Nakajima, A. Tanaka, Y. 

Takano, H. B. Wang, M. Nagao, H. Hatano, and S. Kawa- 

bata, Phys. Rev. Lett. 95, 107005 (2005). 

X. Y. Jin, J. Lisenfeld, Y. Koval, A. Lukashenko, A. 

V. Ustinov, and P. Miiller, Phys. Rev. Lett. 96, 177003 

(2006). 

J. M. Martinis, S. Nam, J. Aumentado, and C. Urbina, 
Phys. Rev. Lett. 89, 117901 (2002); R. McDermott et al., 
Science 307, 1299 (2005) 

E. Chow, P. Delsing, and D. B. Haviland, Phys. Rev. 
Lett. 81, 204 (1998). 

R. Kleiner, F. Steinmeyer, G. Kunkel, and P. Miiller, 
Phys. Rev. Lett. 68, 2394 (1992); R. Kleiner and P. 
Miiller, Phys. Rev. B 49, 1327 (1994). 
T. Koyama and M. Tachiki, Phys. Rev. B, 54, 16183 
(1996); A. Gurevich and M. Tachiki, Phys. Rev. Lett., 
83, 183 (1999). 

S. Coleman, Phys. Rev. D 15, 2929 (1977);C. G. Callan, 
Jr. and S. Coleman, Phys. Rev. D 16, 1762 (1977). 
I. Affleck, Phys. Rev. Lett. 46, 388 (1981) 
G-L. Ingold in Quantum Transport and Dissipation, Th. 
Dittrich, G-L. Ingold, G. Schon, P. Hanggi, B. Kramer, 
and W. Zwerger, , Wiley- VCH (1998). 



